Effect of microscopic disorder on magnetic properties of metamaterials 
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We analyze the effect of microscopic disorder on the macroscopic properties of composite meta- 
materials and study how weak statistically independent fluctuations of the parameters of structure 
elements can modify their collective magnetic response and left-handed properties. We demonstrate 
that even a weak microscopic disorder may lead to a substantial modification of the metamaterial 
magnetic properties, and f0% deviation in the parameters of the microscopic resonant elements 
may lead to a substantial suppression of wave propagation in a wide frequency range. A noticeable 
suppression occurs also if more than 10% of resonant magnetic elements possess strongly different 
properties, but in this latter case the defects can create an additional weak resonant line. These re- 
sults are of the key importance for characterizing and optimizing the novel composite metamaterials 
with left-handed properties at terahertz and optical frequencies. 
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I. INTRODUCTION 

Fabricated composite conductive structures for electro- 
magnetic waves or metamaterials acquire a growing at- 
tention of researches during the last years due to their 
unique properties of magnetic permeability and left- 
handed wave propagation. Being composed of three- 
dimensional arrays of identical conductive elements, the 
metamaterials have much in common with conventional 
optical crystals scaled to support the propagation of mi- 
crowave or Terahertz radiation. In contrast to crystals, 
metamaterials allow tailoring their macroscopic proper- 
ties by adjusting the type and geometry of the structural 
elements. In particular, it is possible to obtain negative 
permeability values in magnetically resonant metamate- 
rials at the frequencies up to hundreds of Terahertz PJQ- 
This appears to be especially useful for the practical re- 
alization of the negative refraction phenomenon Q . 

The simplest method to create a magnetically resonant 
metamaterial is to assemble a periodic mesh of the res- 
onant conductive elements (RCEs), where each element 
is much smaller than the radiation wavelength, and it 
can be well approximated by a linear LC contour. In- 
troduction of a small slit provides the contour with a 
certain capacitance while its shape determines the self- 
inductance. As a result, a resonance of induced currents 
and corresponding magnetization resonance occurs. Re- 
markably, the resistive losses in RCEs are low enough to 
provide the quality factor of the resonance of the order 
of 10 3 0. 

During recent years, a number of ideas have been sug- 
ested for achieving better characteristics Q, tunability 
IE 1 nonlinear wave coupling 0, 0, in composite 
metamaterials by inserting different types of active and 
passive electronic complements into the resonant circuits. 
One of the common features of these seemingly different 
approaches is to modify the macroscopic properties by 
identical insertions into the microscopic structure of the 



metamaterial assuming a technologically ideal fabrica- 
tion process. It was always assumed that all RCEs are 
identical, and no analysis of the effect of random devi- 
ations or disorder was carried out. However, we should 
expect that even weak parameter fluctuations in the mi- 
croscopic parameters may become critical near the mag- 
netic resonance. On the other hand, after a certain op- 
eration time, a small amount of RCEs could experience 
a breakdown operating in a way different from other ele- 
ments. Thus, the problem of disorder appears naturally 
in the theory of composite metamaterials. 

Recently, the study of defect elements on transmission 
properties was studied in experiments with one- and two- 
dimensional metamaterial structures The very first 
study of the effect of disorder in the microscopic structure 
of metamaterials has been made in Ref. [l2|> where the 
magnetic susceptibility for a spatially uniform system fl3l | 
has been averaged with respect to random variations of 
the RCE resonant frequency, and the resulting change 
of the frequency dispersion of the left-handed composite 
system has been found analytically. The method used 
in Ref. |l2j is based on the macroscopic averaging per- 
formed prior to the statistical averaging. This is possible 
only under the assumption that the resonant frequency 
is a slowly varying random function of coordinate and its 
correlation length r c satisfies the inequality 



<C r c < A, 



(1) 



where r resp is the characteristic length of the local re- 
sponse (which is usually of the order of a few lattice con- 
stants 01), and A is the wavelength of the electromag- 
netic radiation propagating in the metamaterial. How- 
ever, the opposite case 



< r 



rcsp i 
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appears to be more realistic because even the nearest 
neighboring RCEs are statistically independent. Then 
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the primary characteristic of the problem is not the aver- 
age susceptibility itself but the current distribution, and 
the magnetic properties of the disordered system are de- 
termined by a macroscopic average of the current. If 
the averaging length is large with respect to the correla- 
tion length, i.e., if the inequality J5J is fulfilled, then the 
macroscopic current is a self-averaged quantity [T^j , and 
it should coincide with its ensemble average. In this case, 
the magnetic properties of the system are described by a 
statistical mean current. 

Below we study systematically the effect of disorder 
on the averaged characteristics and susceptibility of the 
metamatcrials, and consider two practically important 
models of the disordered composite metamaterials assum- 
ing the capacitances C n of different RCEs to be random 
quantities. In the first model, we assume that the capaci- 
tances are completely uncorrelated, i.e. the inequality @ 
is automatically satisfied, but the fluctuations are weak. 
The second model corresponds to small volume density 
n of inserted defect RCEs acting as impurities with sub- 
stantially different capacitance C. The difference can be 
very strong here covering two practically important cases 
of casual RCE breakdown, C — ► oo, and absence of some 
RCEs, C — > 0. The impurities make the medium mi- 
croscopically inhomogeneous. The volume in which the 
local response is formed has to contain a large number 
of impurities. In accordance with the condition ©, this 
implies an additional but not very restrictive condition: 
concentration volume density of impurities should not be 
extremely small, n ^> A~ 3 . 

The paper is organized as follows. In Sec.HHwe develop 
an appropriate theoretical method based on the response 
function of the discrete composite medium , which allows 
describing the properties and response of rather general 
microscopically inhomogeneous media. Section IlIII is de- 
voted to the study of effects of small deviations in each 
resonant element. In Sec. IIVI we deal with the case of 
strong but rarefied impurities. Our general conclusions 
are accompanied by some particular examples calculated 
numerically for the typical metamaterial parameters. Fi- 
nally, in Sec. we discuss the quality and reliability re- 
quirements for electronic components to be used as RCE 
insertions in the composite magnetic metamatcrials. 



II. RESPONSE FUNCTION OF MAGNETIC 
METAMATERIALS 



First we consider an ideal composite metamaterial cre- 
ated by a three-dimensional lattice of identical RCEs. 
The RCEs are placed in the parallel planes normal to 
the z axis (see Fig. 1), and they form a three-dimensional 
structure. We denote N the total macroscopic number of 
the lattice sites with the index n. For the external fields 
oscillating with the frequency w, the relation between the 
external electromotive forces, £ n , applied to RCEs, and 
the induced currents / n is given by the mutual impedance 
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FIG. 1: Schematic of two layers of the metamaterial with ran- 
dom insertions of defect resonant elements (shown in green). 



matrix Il5l Z with N 2 elements: 



^ ^ Znm Im 



£n- 



(3) 



where Z nm = Z(n — m) for a periodic structure. The 
diagonal elements of the matrix Z coincide with the RCE 
sclf-impedance: 



Z(0) = -iujL+ — +R, 

LUC 



(4) 



while the non-diagonal ones are determined by the mu- 
tual inductance: 



Z(n — m) = -iwM(n — m) 



(5) 



We define the discrete Green's function of the com- 
posite metamaterial as the distribution of currents G(n) 
when only one RCE with the index n = is exposed to 
the action of the unitary external electro-motive force, 
i.e. £ n = 5 nt o. Then, Eq. © yields 



^Z(n-m)G(m-l) = 



(0) 



for arbitrary indices n and 1. By the other words, discrete 
Green's function is the inverse matrix of Z: ZG = I. 

The inversion of the difference matrix can be easily per- 
formed via the Fourier transform in the reciprocal space. 
The latter consists of N wave vectors k within the first 
Brillouin zone. We define the Fourier transform /(k) of 
a discrete function /(n) as 



/» = ^E/( n ) e ~ Ikn 



so that the inverse transform has the form 



/(n)=£/(k)e 



For the Fourier spectra, using convolution theorem, 
Eq. © reduces to a simple form Z(k)G(k) = 1, that al- 
lows us to write the discrete Green's function as follows, 
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$>(m)e 



(7) 



The poles of the Green's function [i.e. zeros of the de- 
nominator in Eq. Q] determine the spectrum of linear 
waves that can be excited in the composite mctamaterial. 
In the magnetostatic approximation, only a part of the 
spectrum can be revealed. Conventional relativistic light- 
like electromagnetic waves remain beyond the validity of 
this approach. However, as we demonstrate below, small 
wave vectors and high group velocities of the light-like 
waves make their contribution negligible. The excita- 
tions which determine the Green's function on a micro- 
scopic scale can be well explored in the magnetostatic 
approximation. Being first mentioned in Ref. fl3l | as 
Biot-Savart excitons, these linear magnetostatic excita- 
tions were named soon magnetoinductive waves [Trl Il7j . 
However, these waves can be treated, at least formally, 
as a classical analogue of linear excitations of spin sys- 
tems or magnons, and we use the term magnons below 
in this paper. The spectra of magnons were thoroughly 
studied in one- and two-dimensional structures. The 
spectrum of a three-dimensional structure was explored 
only in the nearest and next-nearest neighbor approxima- 
tions |l6l ITsj , which provide only a qualitative informa- 
tion. As it will be discussed in Sec. IIII 51 below, in order 
to obtain a quantitatively reliable results we should take 
into account hundreds of neighbors. 

After changing the summation over the macroscopic 
number of k- vectors by the integration over the first Bril- 
louin zone, B\, we obtain 



G(n) 



d 3 k 



(27r) 3 n J ^Z(m)t 



-ikm 



(8) 



where n is the volume density of RCEs. With the help 
of Eqs. (JSJ we can rewrite the denominator in the 
following way 



G(n) 



iC 



(2tt) 



n 2 (k) 



d 3 k 



'n J oj* 



ft 2 (k) +iw 7 (k)' 



(9) 



Here the dispersion f2(k) is defined as 



n 2 (k) 



n 

and the resonance width is 



(10) 



7 (k) 



n 2 (k) 

Qua 



(11) 



where ujo — (LC) x l 2 and Q = ujqL/R are the RCE 
resonant frequency and quality factor, respectively. 



To isolate the contribution of singular poles, we split 
the integral in Eq. @ into the real and imaginary parts, 



G(n) 



(27r) 3 n 



[A{p) + iB{p)\ 



(12) 



where the imaginary part allows a straightforward nu- 
merical integrating, 



B(n) 



Q 2 (k) w 2 -Q 2 (k) e' kn , , 
— j-i-^ d 3 k, 13 

[u 2 -n 2 (k)f + Lu 2 j 2 (k) 



whereas the real part 



A(n) 



w ft 2 (k) 7 (k) e 



zkn 



n 2 (k)} A + w 2 7 2 (k) 



d A k 



(14) 



has the contribution at the surface 5k built by 
the magnon wave vectors K(w) obeying the relation 
f2(K(o;)) = u. Since the RCE quality factor is high, 
the singular part dominates in Eq. (|14J) . In the vicinity 
of the surface Sk, we write O(k) ~ to + k ■ v(K), where 
v(K) stands for the group velocity. Next, we integrate 
along the surface normal reducing the volume integration 
in Eq. (|14f> to the surface integral over Sk, 



A(n) 



d 2 K 



S K 



v(K) 



(15) 



which is now more suitable for numerical calculations. 

Analyzing Eqs. (|13fl . 1)15(1 . it is easy to conclude that 
accounting for the light-like modes would not lead to any 
noticeable corrections. These modes are located in the 
center of the Brillouin zone near the point k = and 
their contribution to the integral B is apparently small. 
On the other hand, the light group velocity is about two 
orders of magnitude higher than that of the magnons. 
Therefore, the light-like part in A is also negligible. 



III. WEAK FLUCTUATIONS 

A. Modified magnetic permeability 

In order to model the effect of weak disorder in the 
composite media, we assume that the values of the RCE 
self-impedances experience random uncorrelated devia- 
tions due to the capacitance fluctuations. This may cor- 
respond to the results of a real fabrication process when 
the capacitances of different resonators are not identical, 
and they can be treated as independent quantities. In 
this case the correlation radius coincides with the lat- 
tice constant, r c = a, and the inequality is satis- 
fied. Obviously, strong uncertainty should totally destroy 
the macroscopic mctamaterial response. We assume here 
that the fluctuations are weak and study their effect on 
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the permeability in the second order of the perturbation 
theory with respect to capacitance fluctuations, using the 
methods known from the solid state physics [19 ]. 

We define the local fluctuations of the RCE impedance 

as 



where 



1 

u 



1 

C 



1 1 

a ~ c 



(16) 



(17) 



where angle brackets stand for statistical averaging. To 
calculate the magnetic permeability, we follow the pro- 
cedure described in Ref. 0| an d study the metamate- 
rial exposed to a homogeneous external field H(°> . Then 
Eq. J2J) takes the form, 



^ Z(n - m)J m + i5 n I n = S, 



(18) 



where £ = iiiqujSHz , and S is the RCE area. Iterating 
Eq. ijT%|) . we obtain 

In = h - i ^ G(n - m.)S m I - 

m 

J2G(n-m.)G(m-p)6 m 5 p I p , (19) 



involves the square of the standard deviation 8 2 — (8^). 

After obtaining the effective impedance, we can per- 
form the rest of the procedure described in Ref. 0| and 
finally write the magnetic permeability as follows, 



iZ cS - (2/3)uj(i nS 2 
iZ e s+ (1 / '3)w [IquS 2 



(23) 



We note that the fluctuations contribute to both imag- 
inary and real parts of Z c s- The imaginary part of the 
correction is determined by B(0) and it affects mainly the 
real part of the permeability. The real part of the cor- 
rection leads to an increase of Im(/x zz ), i.e., enhances the 
effective dissipation. Remarkably, this occurs even when 
RCE losses are absent, and the energy is dissipated with- 
out conversion into heating, i.e. this dissipation mecha- 
nism is analogous of the Landau damping. The nature 
of this non-heating dissipation becomes clear if we note 
that it is determined by the term A(0), which combines 
the contributions from all magnetostatic waves excited 
at the given frequency lo. This suggests that the incident 
radiation experiences scattering on the microscopic fluc- 
tuations. Similar effect arising in electromagnetic media 
without positional order (see, e.g. Ref. [2(| and references 
therein) is known as scattering loss or Raleigh scattering. 
In our case, very small correlation radius suppresses the 
scattering into light-like modes, but the scattering into 
magnons is strong. 



where Iq is the current induced in RCEs without param- 
eter fluctuations, 



/„ = £^G(n)4, 



Z 



(20) 



The assumption of weak fluctuations allows us to substi- 
tute Iq instead of I p into the last term of Eq. (|19fl ■ 

The important step in our subsequent analysis is 
macroscopic averaging of this equation. The size of the 
(macroscopic) volume of averaging should be small with 
respect to the wavelength of the external field. On the 
other hand, this size is much larger than the radius of the 
capacitance correlations. Therefore, from the statistical 
point of view, the volume of averaging can be consid- 
ered infinite. As a result, instead of volume averaging 
the statistical averaging can be performed 0|. Taking 
into account that (<5 n ) = 0, and the capacitances of dif- 
ferent RCEs are statistically independent, (<5 m <5 p ) oc S mp 
we obtain that in the second order of the perturbation 
theory the average current induced in RCEs can be pre- 
sented as [cf. Eq. PU))] 



(In) 



£ 

Zctt 



where the effective impedance 

Z cfi = Z + 6 2 G(0), 



(21) 



(22) 



B. Weak disorder in a typical metamaterial 

As a particular example of our theory, we calculate 
the averaged magnetic permittivity numerically for typ- 
ical metamaterial parameters with a weak disorder. We 
assume circular RCEs with radius r$ = 2mm, wire thick- 
ness Z = 0.1 mm, which leads to the self-inductance 
L = 8.36 nHn (see Ref. 0). To obtain RCEs with the 
resonant frequency loq — 6ir ■ 10 9 rad/s (vq = 3GHz), wc 
take C — 0.34 pF. The lattice constants are a = 2.1ro in 
the plane and b = in the z direction. The RCE qual- 
ity factor Q can reach the values of 10 3 pj. However, we 
expect that insertion of diodes or other electronic com- 
ponents can lower this value to Q = 10 2 . 

First, we calculate the linear spectrum of magnon 
waves from Eq. (|10fl and find a strong evidence of long- 
range interaction effects: the lattice sum converges rather 
slow. In particular, in order to get an accuracy of a few 
percent, we have to expand the summation radius to at 
least ten lattice constants. Further increase of the sum- 
mation limits would be unjustified since the sum should 
be calculated over the distances much smaller than the 
wavelength of radiation. We believe that the problem 
of exact calculation of magnon spectrum in the three- 
dimensional case needs a separate detailed consideration. 
On the other hand, we observe that the maximum un- 
certainty in the spectrum takes place along specific di- 
rections perpendicular to the edges of the Brillouin zone. 
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FIG. 2: Real (a) and imaginary (b) parts of the magnetic 
permeability of metamaterial without fluctuations and with 
weak capacitance fluctuations, 5loqC, marked on the plots. 



Therefore, the resulting error in the integrals determining 
the discrete Green's function is extremely small. 

Evaluating numerically the integrals A(0) and B(0) 
according to Eas. l|13|) . I|15|l we obtain the effective 
impedance and corresponding magnetic permeability. 
The frequency dispersion of the permeability fi zz is pre- 
sented in Fig. [21 for the unperturbed metamaterial and 
for several values of the standard deviation. Apparently, 
already a 10% uncertainty changes dramatically the per- 
meability frequency dispersion near the resonance. The 
effect is most pronounced at the frequencies below the 
resonance. In a wide range, the imaginary part of /.t be- 
comes comparable with the real part, and the weakly 
disordered nondissipative medium becomes strongly dis- 
sipative. In the range of negative \i above the resonance, 
the losses are also considerably higher than in a perfect 
metamaterial. As a result, the frequency range appropri- 
ate for the negative refraction shrinks. 



IV. RAREFIED STRONG DEFECTS 
A. Concentration expansion 

We consider now the metamaterial with a small 
amount of randomly distributed defective RCEs which 
differ strongly from the regular RCEs, and therefore can 
be treated as impurities. We assume that the dimen- 
sionless concentration of impurities c = n/n is low, i.e. 
c <C 1. As above, we assume that the deviations from 
the structure parameters appear due to the difference 
in RCE capacitances, and the capacitance of the impu- 
rity, C, can even become very large corresponding to a 
casual breakdown. We neglect any correlation between 
the impurities but assume that two impurities cannot 
be placed on the same site. The applicability condition, 
?i ^> A~ 3 , defines the lower limit for the impurity concen- 
tration, c m in ~ 10~ 3 for a typical metamaterial. Above 
this limit, we can calculate statistically averaged current 
and construct its concentration expansion using standard 
techniques [T3 . l2l | . 

In the system with impurities, the microscopic current 
distribution is substantially inhomogeneous. We are in- 
terested in the normalized averaged value defined as 



(I) 



(24) 



where the summation is performed over a volume con- 
taining a macroscopic number of impurities. The concen- 
tration expansion of the averaged current can be written 
in the following form l2l"| 



(7}=7o + c^[7i( Pl )-7o] 



Pi 



c 2 

y [ j 2(pi,P2)-7i(pi)-7i(P2) + 7o] 



Pl#P2 

Here 



MpO^E 1 ^) 



(25) 



(26) 



is the averaged current if a single impurity is located at 
the site pi, while 

ia(pi,P2) = ^^4 2) (Pi,P 2 ) 



is that for two impurities placed at the sites Pi, P2- The 
terms of higher order involve the averaged solutions for 
more impurities. Clearly, we can write formally the ex- 
pressions for any term of the expansion plj . However, 
in this paper we focus on the strongest effects which are 
linear in the impurity concentration. To calculate the 
corresponding coefficient, we should find first the corre- 
sponding current distribution. 
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The impedance matrix equation |yj in an impure sys- 
tem exposed to a homogeneous external field takes the 
form 



Y J Z(n-m)I m + iA^2s n , p I P = E, 



(27) 



where the summation in the second term in the l.h.s. is 
performed over the impurity sites p only, and 



w \C C 



(28) 



Inverting the matrix Z, we rewrite Eq. I|27|) as follows, 

/ n = / -iA^G(n-p)/ p . (29) 
p 

If a single impurity is located at the site pi, Eq. I|29|) 
taken at n = pi yields 



/n 



Io 



Pl 1 + «AG(0) 
and the current distribution in this case is 



4 1} (PiWo 



f 



iAG(n - pi) 
l-MAG(O) 



(30) 



(31) 



Averaging this expression according to Eq. (|26|l and sub- 
stituting the result into Eq. I|25(l leads to the same results 
(I2I|I and IL'-il) for the averaged current and the mag- 
netic permeability correspondingly, where the effective 
impedance is 



Z e s — Z + c 



iA 



l-HAG(O)' 



(32) 
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FIG. 3: (a) Real and (b) imaginary parts of the magnetic 
permeability of metamaterials without and with impurities 
(impurity concentration is shown near the curves). 



B. Effective magnetic permeability 

To illustrate our results with a specific example, we 
take typical parameters of metamaterials used above in 
Sec. IIII Bl First, we study the effect of impurities with 
an infinite capacitance, C — > oo, which corresponds, for 
instance, to a casual breakdown of varactor diode inser- 
tions. Next, we deal with other impurities which cither 
do not contribute to the magnetization or are just absent; 
this latter case is modelled by setting C — > 0. Although 
these cases are two opposite extreme limits, the result- 
ing averaged permeabilities look similar and differ only 
by few percents. This result is not surprising, because in 
both cases the currents induced in the defect RCEs are 
either small or vanish. 

In Fig. [3] we present the real and imaginary parts of 
the magnetic permeability for several concentrations of 
missing RCEs. As follows from those results, the res- 
onance is attenuated much more than one could expect 
from the fact that a small part of RCEs do not contribute 
to the resonance. The reasons for that enhancement are 



scattering losses which decrease the effective quality fac- 
tor by an order of magnitude or more. We see also that 
the metamaterial is more sensitive to such perturbations 
in the frequency range of negative [l. The composite 
medium can tolerate 5% or less malfunctioning RCEs, 
but for 10% and more a substantial damping of the reso- 
nance occurs, and the imaginary part of the permeability 
becomes comparable with the real one. 

The impurities with a finite resonant frequency, loq ^ 
0, which strongly differs from the main frequency log, 
cause observable effects even if their concentration is ex- 
tremely low. As seen in Fig. 01 only 1% of such im- 
purities can build up their own weak resonance, with the 
position of the resonance shifted from the impurity eigen- 
frcquency. 

In the linear approximation used here, the effect of mu- 
tual interaction of impurities is neglected. Therefore, we 
are not accounting for the magnons arising at frequencies 
close to the impurity resonance. We expect that the cor- 
responding scattering losses will broaden the resonance. 
However, the detailed analysis of this and related phc- 
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FIG. 4: (color online) Real (solid) and imaginary (dashed) 
parts of the permeability of metamaterial with 1% of ran- 
domly distributed defects with the eigenfrequency Cjq — 0.5u>o- 
The inset shows a blow up of the impurity resonance. 



nomcna is beyond the scope of this paper. 



V. CONCLUDING REMARKS 

We have demonstrated that even a weak microscopic 
disorder in the conducting elements of composite meta- 
materials can lead to a substantial modification of their 
resonant magnetic response. According to our results, 
already 10% disorder in the parameters causes substan- 
tial scattering of the incident radiation into magnetiza- 
tion waves. This modifies strongly the macroscopic per- 
meability of the metamaterial leading to the increase of 
losses. We believe the study of disorder and effective 
losses is of a critical importance for engineering novel 
metamaterials with various electronic elements inserted 
into their microscopic resonators. The inserted elements 
should possess no more than a few percent deviation of 



their capacitance and/or inductance to ensure that the 
metamaterial properties are not distorted dramatically. 
Another restriction concerns the insertion stability. The 
data obtained demonstrate that the metamaterial toler- 
ates no more than 10% disorder in the parameters of the 
resonant conductive elements. Casual breakdown of a 
greater part of the insertions causes strong damping of 
the wave propagation in metamaterials. 

Our results that even a small amount of defects can 
build up a noticeable additional magnetic resonance look 
useful for suggesting a simple methods for the sensitive 
quality control of metamaterials. If an RCE contains not 
a single capacitive element but a combination of inser- 
tions with an ordinal slit, a casual insertion breakdown 
would switch the RCE to another resonant frequency. 
The corresponding narrow gap (or peak) in the metama- 
terial transmission can be used as a sensitive indicator of 
the concentration of damaged RCEs. 

Finally, we point out that the effects of microscopic 
disorder would be crucially important in the nanostruc- 
tured metamaterials even with simple RCEs without elec- 
tronic components. Clearly, the fabrication of resonant 
elements on such scales is less accurate, so that random 
fluctuations of the RCE shapes can lead to deviations of 
the self-inductance and capacitance. Disorder in posi- 
tion of RCEs determines the deviations of mutual induc- 
tance. We expect that the results obtained here can be 
important for qualitative estimations of nanostructured 
left-handed media. However, the problem requires a sep- 
arate and more systematic analysis. 
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